Ramping fermions in optical lattices across a Feshbach resonance 
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We study the properties of ultracold Fermi gases in a three-dimensional optical lattice when 
crossing a Feshbach resonance. By using a zero-temperature formalism, we show that three-body 
processes are enhanced in a lattice system in comparison to the continuum case. This poses one pos- 
sible explanation for the short molecule lifetimes found when decreasing the magnetic field across a 
Feshbach resonance. Effects of finite temperatures on the molecule formation rates are also discussed 
by computing the fraction of double-occupied sites. Our results show that current experiments are 
performed at temperatures considerably higher than expected: lower temperatures are required for 
fermionic systems to be used to simulate quantum Hamiltonians. In addition, by relating the double 
occupancy of the lattice to the temperature, we provide a means for thermometry in fermionic lattice 
systems, previously not accessible experimentally. The effects of ramping a filled lowest band across 
a Feshbach resonance when increasing the magnetic field are also discussed: fermions are lifted into 
higher bands due to entanglement of Bloch states, in good agreement with recent experiments. 

PACS numbers: 03.75.Ss, 05.30.Fk, 71.10. Ca 
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I. INTRODUCTION 

Ultracold atoms loaded into optical lattices are nearly 
ideal experimental realizations of nonrelativistic quan- 
tum lattice models (Jj . Experiments on bosonic atoms in 
optical lattices have demonstrated the transition from su- 
perfluid to Mott insulator states , and probed momen- 
tum distribution and excitation spectra in strongly inter- 
acting bosonic systems 13, Q . Comparisons to quantum 
Monte Carlo 0, IE IS an d density matrix renormal- 
ization group calculations exhibit good quantitative 
agreement and confirm the quantitative applicability of 
simple quantum lattice models to describe atomic gases 
in optical lattices. Bosons in optical lattices thus provide 
an interesting test bed for theoretical models, such as the 
recently introduced ring-exchange model [ill). 

In contrast to bosonic systems, the physics of strongly 
interacting fermions in two- and three-dimensional lat- 
tices is not yet fully understood. Indeed, accurate nu- 
merical simulations on interacting fermionic systems in 
two- and three-dimensional lattice models are nondeter- 
ministic polynomial hard |l2^ . and no general numerical 
or analytical solution exists. Experiments on ultracold 
fermionic gases in optical lattices could thus be very use- 
ful in elucidating the properties of the fermionic Hubbard 
model [T3lll4j and to probe exotic quantum phases, such 
as d-wave resonating valence bond phases [l5| . 

Progress in cooling techniques for fermionic atoms has 
recently allowed the first experiments on cold fermionic 
gases in optical lattices to be performed by T. Esslinger 
and coworkers [Tfjl IrH ITsI ]. In a first experiment they 
loaded a weakly interacting cold Fermi gas into an op- 
tical lattice and observed the momentum distribution 
function. Ramping this Fermi gas across a Feshbach res- 
onance by increasing the magnetic field has been seen to 
lift a fraction of the atoms into higher bands. In a second 
experiment, the atoms were ramped across the Feshbach 
resonance by decreasing the magnetic field, which leads 



to the formation of molecules jl9| . 

In this paper we present analytical and numerical cal- 
culations on fermions in three-dimensional optical lat- 
tices to quantitatively explain the experiments of Kohl 
et al. 01 when ramping up the field. We start in SecHTI 
with a calculation of the two- and three-particle densities 
in optical lattices to compare loss rates in the continuum 
to lattices. Next we consider the effects of finite temper- 
atures in Sec. Illll and relate the double occupancy in the 
lattice to the temperature of the system, thus providing 
a means for thermometry |2fj| in fermionic experiments. 
Our results show that at experimentally accessible tem- 
peratures, the double occupancy of lattice sites is sub- 
stantially reduced compared to the ground state; this has 
important consequences for the quantitative explanation 
of the experiments of Ref. 0], discussed in Sec. IIVI 



II. TWO- AND THREE-PARTICLE DENSITIES 

We start by solving the Schrodinger equation for non- 
interacting fermions in an optical lattice, in order to com- 
pute the two- and three-particle densities. 



A. Schrodinger equation 

Superimposing two counterpropagating running-wave 
laser beams of wavelength A = 27r/fc (k is the wave vec- 
tor) in all three spatial directions yields a static potential 
U(r) with periodicity a = A/2: 



U(r) =^[/sin 2 (fc l x l 



(1) 



i=l 



where r = (xx,X2,X3). We consider a cubic system 
with iV 3 sites, and a linear extent L = aN. In or- 
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FIG. 1: (Color online) Left: Five lowest bands in the first 
Brillouin zone plotted as a function of the wave vector p = 
(p, 0,0) for V = 10 (solid lines) and 40 (dashed lines). Right: 
Band gap Ao,i between the lowest and first excited bands as 
a function of the potential V in units of E-r. 



der to describe fermions in an optical lattice, the time- 
independent Schrodinger equation 



<F(r) = EV(r) 



(2) 



has to be solved for each spin state. Here m repre- 
sents the atomic mass. Using a product ansatz \l/(r) = 
i'i{xi)4'2{x2)^z{xz) with E = Ei + E 2 + E 3 , each com- 
ponent can be solved via Bloch states 



with periodic functions 



a) 



(Xi) 



(3) 



(4) 



and corresponding energy eigenvalues E ni jPi . Here rij € 
N is the band index and p, e I the quasimomentum. 
In order to describe finite lattices, we impose periodic 
boundary conditions [ipi,m,Pi( x i) — ^i,ns,p 4 { x i + L)] to 
the solution, yielding a discrete set of N allowed quasi- 
momenta pi = —ir/a, —ir/a + 2tt / L . . . , ir/a — 2?t/L in 
the first Brillouin zone (first BZ) given by the interval 
(— 7r/a, 7t/a). A full three-dimensional Bloch state ^ n ,p 
is thus parametrized by two vectors n = (n-i, rc.2, n.3) 
and p = (pi,P2iPs)i a s well as its energy E np — 
Em,pi + En 2 ,p 2 + ^"3,p3i where each band carries iV 3 
allowed states for each spin state. Figure ^ left panel, 
shows the five lowest bands for V = U /E-r = 10 and 40, 
where Er = h 2 k 2 /2m. The right panel of Fig. ^ shows 
the band gap Ao^ between the ground state and first 
excited state as a function of the potential V. Inserting 
Eq. iJSJl into Eq. ©, we obtain 



1 

k 2 



dx 2 



2ip 



d_ 

dx 



Vsin (kx) 
= e 



U n ,p{x) (5) 
n,p'U"n,p{x) . 



In Eq. © we have expressed the potential U in units of 
the recoil energy En = h 2 k 2 /2m, i.e., V = U/En and 



e UjP = E ny p/E-n. We solve Eq. JSJ with a Fourier ansatz 
for periodic functions 



+c 



i lP (x) = lim Un >Pl j 



j=-c 



a/2 



-a/2 



(6) 



(7) 



Note that in Eqs. JSJ and JJJ we introduce a cutoff c. If 
we consider only the lowest five bands e n ,p, n = 1, ... 5, 
convergence beyond machine precision is attained already 
for c = 6. 



B. Two- and three-particle densities 

Experimentally, molecule formation rates [Tfij ] have not 
been as high as expected. This can be attributed either 
to temperature fluctuations or to three-body processes 
hampering the molecule formation. Naively one would 
expect that increasing the potential depth V would sup- 
press three-body processes disrupting molecule forma- 
tion. In order to test this, we compute the three particle 
density and show that at moderate potential depths V, 
in a range relevant for experiments, it is actually larger 
than without an optical lattice. 

To calculate the probability of having three fermions 
in a small volume, we start with a filled lowest band 
n = no = (0,0,0), i.e., each site carries two particles 
with opposite spin states. The ground state |4> ) is then 
given by 



= n «i,p. a,p,i i°> 

pGfirstBZ 



(8) 



where s (a niP . s ) creates (annihilates) a particle in a 
Bloch state l^n.p) with spin s E {T, !}• The product 
over p is carried out over the allowed states in the first 
Brillouin zone. The corresponding field operators defined 
in terms of Bloch states are given by 



p,s > 



$ s (r) = ^*„, p (r)a n , 

n,p 

$t( r ) = £** >p (r)at, p>s , 



(9) 



n.p 



where the summation is carried out over all band indices 
n e Nq and allowed p in the first Brillouin zone. The 
three-particle density ps in a small volume at a given 
lattice site is then given by 

P3 (r, r') = ($ 1 $} (r)${ (r)*{ (r')* T (r')*l (r)* T ( r ) I $0) , 

(10) 

where we have one up and one down spin at position 
r, and one up spin at position r'. The operators W|(r) 



3 



[^(r)] create (annihilate) a particle with spin state s £ Eq. (|10|l . After a straightforward derivation one obtains 
{T,l} at a space point r. These can be inserted into for the three-particle density 



ps(r,v') = [l*n , Pl (r)| 2 |* no . P2 (r)| 2 |* no , P3 (r)| 2 - ^ , Pl (r)* no , Pl (r')|* n ^ P2 (r)| 2 ^ o!p3 (r')* no ,p 3 (r)] . (11) 

Pl,P2,P3 



Flipping the signs of the spins at r and r' in Eq. Hl()(l 
does not affect the expression in Eq. (|llfl . 

The two-particle density p 2 for two atoms in the same 
spin state is given by 

p 2 (r,r') = ($ |${(r')$}(r)$ T (r)$ T (r')|$o), (12) 
and can be evaluated in a similar way as Eq. (|l(Jfl : 
p 2 (r,r') = J2 [l*n , Pl (r')| 2 |*„ r , P2 (r)| 2 - (13) 

Pl>P2 

*: 0>P1 ( r ')*n , Pl (r)*; , P2 (r)* no , P2 (r')] . 

Finally we compute the probability p?, of finding two 
fermions in opposite spin states at positions r and r' 

^(r,r') = (*o|*|(r / )*I(r)*i(r)$ T (r , )|*o) (14) 
which can also be expressed as 

p° 2 (r,r') = J2 [l*n„, Pl (r')| 2 |*n , P2 (r)| 2 ] . (15) 

P1.P2 

Because of the short-range nature of the three-body in- 
teractions, we are most interested in the limit |r— r'| — * 0. 
Since both p2 and p^ vanish in that limit due to the Pauli 
principle, we calculate their average over a small volume 
51r (a/2 < R < a/32) around the origin 

~p {v ^w~\l p(°' r V- ( 16 ) 

|"r| Jn R 

In Fig. we show data for the averaged two- and 
three-particle densities according to Eq. (jpp|l as a func- 
tion of the potential depth V. The data show a marked 
peak that exceeds unity for moderate potential depths. 
This leads to three-particle processes that in turn lead 
to molecular loss in the fermion gas. For large potential 
depths both densities are exponentially suppressed. This 
is in contrast to a configuration with only two fermions 
having different spins, as can be seen by divergence of ~p\ 
in Fig. 

C. Partially filled lowest bands at zero temperature 

So far we have not taken into account the Gaussian- 
shaped intensity profile of the lasers generating the opti- 
cal lattice which yields a superimposed harmonic confine- 
ment |2l| . The inhomogeneity of the trapping potential 
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FIG. 2: (Color online) Averaged three-particle density ~p 3 (V) 
and average two-particle density p 2 (V) via Eq. I|16|l as a func- 
tion of the optical potential depth V calculated on cubes of 
linear extent R = a/32 in lattices with TV = 65. The graphs 
are normalized in such a way that ~p 3 (V = 0) = 1 (free fermion 
gas). The dashed horizontal line marks where the densities 
cross unity. As shown in the inset, for a completely filled low- 
est band and V > 6, the probability of finding two particles 
per site in the same spin state is exponentially suppressed, 
whereas the probability of finding two fermions in different 
spin states, 'p^iV), diverges. 



influences the spatial density distribution of the parti- 
cles in the lattice. The effects of the trapping potential 
are often considered in a "local density" approximation, 
where the local properties in the trap are approximated 
by results obtained from a homogeneous system with the 
same local density. 

The partially filled ground state in terms of Bloch 
states is given by 

= II «Lp A, P a 1°) ■ 

peA 

where A is the region of the first Brillouin zone inside the 
Fermi surface. We can again compute the three-particle 
density now restricted to momenta inside A and obtain 



4 



p 3 (r,r')= Xa(Pi)xa(p 2 )xa(p 3 ) [|*no, Pl (r)| 2 |*no,p 2 (r)| 2 |*n„,p3(r) 



^o,P I ( r )^o,p 1 (rOl*n ,p 2 (r)| 2 ** , P3 (r')* I10 ,p 3 (r)] , (18) 
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FIG. 3: (Color online) Averaged three-particle density p 3 (V) 
via Eq. 1161 as a function of the optical potential depth V 
for R — a/2. The graphs are normalized in such a way 
that p 3 (V = 0) = 1 (free fermion gas). The data are for 
a partially filled lowest band in the ground state. The three- 
particle density decreases monotonically when the filling fac- 
tor pi = iVo/(2iV 3 ) is decreased (data for pi = 0.1, 0.2, 0.3, 
...1.0). This shows that three-particle processes are sup- 
pressed when moving away from the band-insulating regime. 
The inset shows a log-log plot of the averaged three-particle 
density divided by pi, J>^iy)/p\, as a function of the filling 
factor pi for different values of the dimensionless potential 
depth V. For all V, ~p 3 {V)/p\ decays according to a power 
law in pi. The dashed line in the inset is a guide to the eye 
proportional to p?' 391 . The data shown are for a — 1 and 
N = 65. 



where the characteristic function xa is 



/ \ = / 1, P e A, 
Xa [P) - | otherwise. 



(19) 



In Fig. |31 the averaged three-particle density for a par- 
tially filled lowest band is plotted as a function of the 
potential strength V for several filling fractions pi = 
Nq/(2N 3 ), where p\ = 1 for a completely filled lowest 
band. For all values of V, the three-particle density de- 
creases with a power-law behavior in p\ (inset of Fig.|3J), 
thus showing that in parabolic traps three-body processes 
can be suppressed in partially filled lowest bands due to 
the reduced filling of the lattice. 



III. THERMOMETRY AND EFFECTS OF 
FINITE TEMPERATURES 

In order to obtain a realistic estimate of the molecule 
formation rate, as well as the lifting of fermions into 
higher bands (cf. Sec. IIV|) the effects of finite tempera- 
tures must be taken into account. We start from a tight- 
binding-like Hamiltonian for the fermions, 



r 2 n 



(20) 



Cj pCj.a is the number 



where the optical lattice potential is quadratic (r 2 
y 2 + z 2 ) in units of A/2 and rij t 
operator. 

Using experimentally relevant parameters |l6j | for the 
trap frequencies, ui x — 2tt x 211 Hz, u> y = 2ir x 257 
Hz, uj z — 2tt x 257 Hz, which corresponds to a potential 
depth of V[E R ] = 15.64, we obtain t/E R = 0.00652 [13, 
K Z /E R = 0.00411, and K X /E R = K y /E R = 0.00609. 
We find that a system of size 100 3 is large enough to en- 
sure that the particle density is zero at the edges. For a 
given temperature in units of the Fermi energy E-p, we 
tune the chemical potential such that for T = 0, as well 
as for a given T > 0, we obtain 40000 particles in the 
trap [l6j| . 

The fraction d,2 of double-occupied sites is given by 



do = 



Ei<"*T> 



(21) 



Here rij^ is the number operator for particles with spin a. 
The second equality in Eq. I|21|) follows from the fact that 
the particles are noninteracting. Using pj — (rij ia ) for the 
expectation value of the number operator we obtain 

, E, p] 



E, p 3 



(22) 



Data for d,2 as a function of temperature T are shown in 
Fig. 0] for the experimental parameters used in Ref. |16| . 
Empirically we determine for several sets of experimental 
parameters that in the band-insulating regime for low 
temperatures and V > 5 



d 2 ps a + 6exp(— jT), 



(23) 



where the temperature is measured in units of the Fermi 
energy Ep. Note that for larger V we obtain a = and 
b = 1 in Eq. (|23|l . i.e., c?2 ~ exp(— 7T). This can be seen 
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FIG. 4: (Color online) Fraction of double-occupied sites in 
an optical lattice of size 100 3 with 40000 particles for the 
experimental parameters presented in Ref. |h[|. The tem- 
perature is measured in units of the Fermi energy Ey- The 
fraction of double-occupied sites decays exponentially with in- 
creasing temperature. In the experiments of Kohl et al. [Tfij | 
d,2 = 0.43 ± 0.05 (dashed ± dotted lines). This means that 
the experiment is done at T/E-p ~ 0.47 ± 0.07, which is con- 
siderably hotter than expected. The inset shows the fraction 
of double-occupied sites as a function of the potential depth 
at T = 0. 



in the inset to Fig. 0] for deep lattices, i.e., V — > oo, 
we obtain d,2(T = 0) — > 1 (data for 40000 particles and 
the experimental parameters used above). For the ex- 
perimental parameters of Kohl et al, we find 7 = 1.818 
(a = 0, b = 1) for < T/E F < 0.65. Because the double 
occupancy depends only on the Fermi distribution and 
the density of states (which in turn depends on the space 
dimension), we expect the functional form to be indepen- 
dent of the trap parameters and only to depend on the 
space dimension of the lattice (which affects the density 
of states). Thus the molecule formation rate can be used 
in experiments to determine the temperature of the sys- 
tem, something which previously has not been accessible 
experimentally. These, as well as similar results obtained 
experimentally by Stoferle et al. 18] show that current 
experiments with fermions in optical lattices are done at 
temperatures considerably higher than in the continuum 

(Fig.g. 



IV. RAMPING ACROSS A FESHBACH 
RESONANCE 

A magnetic field tuned near a Feshbach resonance 
[23| can be used to vary the two-body s-wave scattering 
length a between two particles in different spin states 
over a broad range which at low temperatures becomes 



the only relevant quantity to describe interactions. Ex- 
perimentally [Til the s-wave interaction is varied on a 
short time scale compared to the tunneling time between 
two adjacent potential wells. Thus, the system in the 
band insulator regime can be regarded as an array of de- 
coupled harmonic potential wells. Therefore, in the pres- 
ence of a magnetic field, a cell of the optical lattice with a 
completely filled lowest band is well described by a two- 
body problem in a harmonic potential with zero-range 
interaction, a problem that has been studied theoreti- 
cally in different approximation schemes l24l I25L |26L Wf\ . 

The scenario considered by Kohl et al. [jjj is to start 
with a completely filled lowest band at a certain exter- 
nal magnetic field B\ below the Feshbach resonance. By 
crossing the resonance with a given velocity (in the ex- 
periment the ramping velocity is determined by maxi- 
mizing the molecule formation rate as a function of the 
ramping velocity) and ending at a final magnetic field B{, 
the large variation of the scattering length ao causes an 
interaction-induced transition between Bloch bands, i.e., 
a fraction of fermions is transferred to higher bands. We 
can understand this problem theoretically by restricting 
it at first to two particles situated at ri and r 2 , and 
choosing Bi and Bf such that ao(-Bi) = ao(Bt) = 0. 
The general solution to this two-body problem [24| is a 
wave function 0(r, R) expressed in ccntcr-of-mass R = 
(ri + r 2 )/\/2 and relative coordinates r = (ri — r 2 )/v / 2- 
For a vanishing scattering length ao, the initial and final 
wave function 0;/f (r, R) can be expressed as 



9 i/f (r,R) = $ ni/f (R)0„ i/f (r) 



(24) 



where $ n (R) is the harmonic oscillator eigenstate with 
quantum number n £ N 3 and 



(j> n (r) = 7T 



= ^-3/4 



LT (0) 



-1/2 



2/2 LT (r 2 ) (25) 



is the harmonic oscillator s-wave function with quantum 
number n given in terms of generalized Laguerre polyno- 
mials l\[ 2 (r) j2^| . The relative part of the wave function 
[cf. Eq. (|24|l ] . i.e., the oscillator s-wave function, is then 
shifted upward by one quantum number when crossing 
the Feshbach resonance: 



m =m + l, 



(26) 



whereas the quantum number of the center-of-mass con- 
tribution n remains unchanged. 

In the following we briefly outline our calculations be- 
fore elaborating on the details. We start by expressing 
the Bloch bands ^n.p around a potential minimum R by 
a product of two localized Wannier states W a< R(r), i.e., 



VF n , R (r) = 



1 



ip-R 



jV 3 / 2 ^ 
P 



^n, P (r), 



(27) 



where the summation is performed over the iV 3 momenta 
p in the first Brillouin zone. While the Bloch states are 
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extended over the lattice, the Wannier functions peak 
at a given lattice site with weight in neighboring sites 
being exponentially suppressed. The rate of exponen- 
tial decay of Wannier functions in one dimension for a 
given energy band is related to the energ y g ap between 
that particular band and any other band [29j. Although 
in one space dimension the energy bands can always be 
well separated by increasing the optical potential depth 
V, the situation in three dimensions is more compli- 
cated due to entangled energy bands. If there is level 
crossing it is not possible to construct localized wave 
functions from Bloch states from one of these bands 
alone 30]. In this case localization of the correspond- 
ing Wannier functions is still possible because the poten- 
tial V(r) allows for the wave functions to separate, i.e., 
W n ,n(r) = w nuRl (x 1 )w n2iR2 (x 2 )w n3tR3 (x 3 ), where the 
single factors w nuRi (xi) can always be localized. Note 
that the calculations using Wannier functions are done 
with a cutoff in the Fourier expansion of c = 6. Con- 
vergence tests show that our results are approximately 
independent of the lattice size for TV > 65, and we per- 
form the following calculations for N = 65. 

In a next step, the Wannier functions are approximated 
with harmonic oscillator eigenstates and are split into 
center of mass and relative parts. Each relative part of 
an element in this approximation is shifted upward by 



a quantum number, i.e., 4> n ( r ) ~ * </>n+i(r). The final 
shifted two-particle wave function after passing the Fes- 
hbach resonance is then transformed back to absolute co- 
ordinates where wave functions in different Bloch bands 
then mix. Applying this procedure to all sites in the 
lattice yields the final excited state of all 27V 3 particles. 

In second quantization, we first express the Bloch 
states in Wannier functions by means of a unitary trans- 
formation [inverse of Eq. (|27() ]: 

R 

The operators and £u n R s in Eq. create a par- 

ticle with spin s in a Bloch state ^n.plr) or Wannier 
state W n n(r), respectively, where R denotes a site in 
the lattice. We next expand the Wannier states in har- 
monic oscillator wave functions, up to order o. Thus the 
particle number operator for a Bloch state v I / n,p(r) with 
spin s given by 

/,,,,,/.,,., = ^3 £ < P l " H i -;,H,''n.H , (29) 
R,R' 

can be approximated as 



/n,p,s/n,P,s 



1 

iV 3 



E 

R,R' 



'»• »" " '< r- a 



r »a n ,R,s 



-n,n« n R 



-n,n+u^n+u,R/ ,s 



(30) 



where a n R (a n ,R,s) creates (annihilates) a particle with 
spin s in the harmonic oscillator state $ n (r — R)- Here 
c n.k = ($k|Wn,o) is the projection of the Wannier state 
W n , onto $ k and u € {(±2, 0,0), (0, ±2,0), (0,0, ±2)}, 
provided that for the summations in Eq. (|30|1 only u are 
considered for which the components of u + n are non- 
negative. Note that, because of the s-wave nature of 
the ground state, only harmonic oscillator wave functions 
with an even total order appear in this expansion. 

To first order in the expansion in harmonic oscilla- 
tor wave functions, only terms containing c nin c nn are 
considered whereas in second order terms containing 
c n ,nC n ,n+u are added. In order to measure the quasimo- 
mentum distribution, one starts with the initial many- 
particle wave function 



i*>>=n 



■^no.R.T^lo.R.I 



|0) 



(31) 



R 



with no = (0, 0, 0), where the single w no R jii^ R ■ can 
be expanded in first or second order (in harmonic oscil- 
lator functions) in analogy to the terms of the sum in 
Eq. (|30|l . For a first- and second-order expansion in har- 



monic oscillator states, the transformation to relative and 
center-of-mass motion coordinates can be expressed in a 
form that transforms back to normal coordinates after 
the increment of the s-wave functions. Starting from a 
ground state we obtain 



i*f) = IP«Ui°>> 



(32) 



r 



where R depends only on the site R and the expansion 
order o € {1,2}. The lowest-order result is 



F-i td — ; — C„„ nr , \ ~ CI — t-> -t CI 



1,R /^no.no I 2 "no ,R,T (0,0,2) ,R, J. 



(33) 



_^2 a (o,o,i),R,T a (o,o,i),R,i + 2 a (°'°' 2 )' R 'T an o,R-4 



1 



at t> 1-al 



1 



2 tx n ,R,T u '(0,2,0),Ra /jj a (0,l,0),R,T a (0,l,0),R,i 

J_-t -t M ~t ' 

v /2" a (l,0,0),R,T a (l,0,0),R,l "r 2 a (2,0,0),R,T a no,R-,l 
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FIG. 5: (Color online) Upper set of curves: Fraction of 
fermions transferred to higher bands for a first- (o = 1) and 
second-order (o = 2) expansion [Eq. 1351 1 plotted as a func- 
tion of the optical potential depth V in units of En- For both 
o 6 {1, 2}, the fraction calculated in this work is nonvanish- 
ing, thus proving an interaction-induced transition between 
lowest bands, in qualitative agreement with results of Kohl 
et al. (16| (dashed line). The lower set of curves takes into 
account that, at finite temperatures, not all sites are double 
occupied. Using the parameters from Kohl et al. we expect 
T « 0.47i?F. The small discrepancy with the experimental 
data can be explained due to the parabolic trapping potential 
(see Sec. lira. 



and the more complex higher-order expression -F^r ^ s 
omitted for the sake of brevity. Note that the quantum 
numbers of the creation and annihilation operators for 
harmonic oscillator states (a"' ) for each term always add 
to an even number because of the s-wave nature of the 
local wave functions. 

The transformed multiparticle wave function Eq. i|32|) 
can now be used for the quasimomentum measurement 



1 



N(o,n,p) ee 5<*f|^ iPtT /^p J T + /:, Pl ;/n,p,il*f>, (34) 

where the particle number operator /i PiS / n ,p,s from 
Eq. (|30|) is also considered in order o. For a first- and 
second-order expansion in terms of harmonic oscillator 
eigenstates, the calculation of the quasimomentum distri- 
bution is carried out analytically and the corresponding 
fraction of particles transferred to higher bands is defined 
as 



?i 



d = 



Epei.'BZ N (°^ P) 

EnEpel.'BZ^(°> n >P) 



(35) 



where the summation over the index n is carried out over 



In Fig. the fraction of particles transferred to higher 
bands from Eq. I|35l) is plotted as a function of the opti- 
cal potential depth V for a calculation in first- (o = 1) 
and second-order (o = 2) expansion. Our results show 
that already for V^-ErJ ~ 15 the data have converged to 
the limit of a deep lattice. Over the whole range of the 
potential depth and for both orders, the fraction is non- 
vanishing, thus proving an interaction-induced transition 
between lowest bands in qualitative agreement with the 
results of Kohl et al. Note that while the relative 

part of the harmonic oscillator wave functions is lifted 
to the first band, this is not the case for the particles in 
the lattice. In the limit of deep lattices we find that [see 
Eq. (23] only 75% of the particles are lifted: 50% into 
the first band and 25% into the second band, while the 
rest remain in the ground state. Finally, in Fig. [5] the 
quasimomentum distribution N(o, n, p) from Eq. (|34|> 
with o = 2 is plotted in an extended zone scheme be- 
fore (left panel) and after (center panel) ramping across 
the Feshbach resonance in comparison to the experimen- 
tal results of Kohl et al. after crossing the resonance 
(right panel). The numerical data are for zero tempera- 
ture and illustrate how higher bands are populated when 
crossing the resonance. In Fig. |7|we show the quasimo- 
mentum distribution in an extended zone scheme plot at 
T/Ef = 0.47 before (left panel) and after (right panel) 
crossing the resonance. The data after crossing the reso- 
nance have been computed by taking into account that at 
T/Ep = 0.47 (th e temperature used in the experiments 
of Kohl et al. only 43% of the sites are double occu- 
pied in the first Brillouin zone. In addition, higher-order 
zones are populated due to thermal fluctuations. The 
theoretical and experimental data agree well. Note that 
Kohl et al. used a trap with direction-dependent frequen- 
cies, which is why the quasimomentum distribution does 
not possess fourfold symmetry. 



V. CONCLUSIONS 

Fermionic systems in optical lattices show rich physical 
phenomena, and the physics of strongly interacting mod- 
els is far from being understood. Before cold fermionic 
gases in optical lattices can be used as quantum simula- 
tors 0, fl3l ITil , the agreement between experimental 
measurements and theoretical results has to be carefully 
checked for solvable models such as the ones discussed 
in this paper. One important issue is to ensure that the 
temperatures at which the experiments are performed 
truly probe the ground-state properties that are desired 
in a quantum simulation. 

The calculations presented here for large three- 
dimensional fermionic systems rely on the fact that in 
experiments |0, 0, Il8f . the initial and final states of 
the fermions are essentially noninteracting. Strong inter- 
actions are present [Tflj only during the short ramping 
periods, and the effect of the interaction could be consid- 
ered locally on the few-body states at each lattice site. 
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FIG. 6: (Color online) Quasimomentum distribution in the (pi,p2) plane, integrated over the third component p$ at T = 0. 
Left: Momentum distribution before ramping across the Feshbach resonance. The first Brillouin zone (black box) is filled. 
Center: Numerical calculation in second order, showing the momentum distribution after crossing the Feshbach resonance. 
Data for a potential depth V[£r] = 40 at which the sets of entangled energy bands are well separated from each other to avoid 
crossovers. The numerical data shown are for a — I and N = 65. Right: Experimental data taken from Kohl et al [l6| . Note 
that in the experiment the optical lattice is turned off before the quasimomentum is measured. 




FIG. 7: (Color online) Quasimomentum distribution in the 
(p\,P2) plane, integrated over the third component pz at 
T ~ 0.47-Ef, the temperature used in the experiments of Kohl 
et al. [lr^ | (i.e., d.2 = 0.43). Left: Momentum distribution 
before ramping across the Feshbach resonance. Right: Nu- 
merical calculation in second order, showing the momentum 
distribution after crossing the Feshbach resonance, taking into 
account thermal effects. 

In order to explain the low molecule formation rates 
observed in experiments, we have computed the two- and 
three-particle densities for fermions in optical lattices. In 
general, one would expect that the three-particle density 
decreases in a lattice system. Our results show that non- 
trivial processes occur when placing a fermionic gas in an 
optical lattice. First, the three-particle density increases 
with increasing potential depth, since up to intermedi- 
ate lattice depths the main effect of the lattice is not 
a localization of the particles on individual lattice sites, 
but rather the reduction of the effective volume acces- 
sible to the fermions. The reduced volume leads to a 
higher fermion density at the lattice sites, and hence to 
larger three-particle densities and scattering. This in- 
crease of three-particle densities, which happens for ex- 
perimentally relevant potential depths, provides an ex- 
planation, besides ramping speeds and imaging quality, 



why the molecule formation rates and lifetimes in exper- 
iments | Hi are lower than expected. The three-particle 
density decreases only for very deep lattices, when the 
Wannier functions become strongly localized on single 
lattice sites and three-body processes should not influ- 
ence the experiments already for F[£e] > 30 (l6j . 

Our calculations show that the fraction of double- 
occupied sites in the optical lattice varies approximately 
exponentially with the temperature. While the direct de- 
termination of the temperature of an experimental setup 
has proven to be difficult as well as inaccurate, the count- 
ing of particles can be done to good precision 0, and 
therefore the fraction of double-occupied sites can be de- 
termined. Our results provide a theoretical foundation to 
recent experimental work by Stoferle et al. [is| in which 
the molecule fraction can be used for thermometry in ex- 
periments, something that was previously not available 
for fermions in optical lattices. In particular, our re- 
sults using standard experimental parameters show that 
experiments are currently performed in a temperature 
regime considerably hotter than for systems in the contin- 
uum. It is furthermore important to perform these mea- 
surements in deep lattices, as otherwise increased three- 
body scattering processes can influence the molecule for- 
mation. 

When tuning the system across a Feshbach resonance 
by increasing the field and thus pushing fermions into 
higher bands we find good quantitative agreement be- 
tween our results and the experiments of Kohl et al. |l6| . 
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